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ABSTRACT 

Aims. Relativistic bending in the vicinity of a massive body is characterized only by the post-Newtonian parameter y within the 
standard parameterized post-Newtonian formalism, which is unity in General Relativity. Aiming at estimating this parameter, we use 
very long baseline interferometry (VLBI) to measure the gravitational deflection of radio waves emitted by distant compact radio 
sources, by Solar System bodies. 

Methods. We analyze geodetic VLBI observations recorded since 1979. We compare estimates of 7 and errors obtained using various 
analysis schemes including global estimations over several time spans and with various Sun elongation cut-off angles, and analysis of 
radio source coordinate time series. 

Results. We arrive at the conclusion that the relativistic parameter 7 cannot be estimated at better than 2 x 10"**. The main factor of 
limitation is the uncertainty in the determination of (global or session-wise) radio source coordinates. A sum of various instrumental 
and modeling errors and analysis strategy defects, that cannot be decorrelated and corrected yet, is at the origin of the limitating noise. 

Key words. Astrometry - Relativity - Techniques: interferometric 
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1. Introduction 

One of the cornerstones of test of General Relativity (GR) is the 
measurement of light deflection in the vicinity of the Sun. In the 
parameterized post-Newtonian (PPN) formalism (Will 1993), 
which contains 10 parameters, the predicted angle of deflection 
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where G is the Newtonian gravitational constant, c the speed 
of light in a vacuum, M the mass of the deflecting body, b the 
impact parameter (defined as the minimal distance of tlie ray 
to tlie center of mass of tlie deflecting body), (f> the elonga- 
tion angle between the deflecting body and the source as viewed 
by the observer and 7 is the PPN parameter characterizing the 
space curvature due to gravity. (See, e.g., Misner et al. 1973, 
Will 1993, and more generally speaking for an axisymmetric 
body, Le Poncin-Lafitte & Teyssandier 2008.) Thus, a grazing 
ray at the Sun's limb is deflected by ~1.7". In GR j - 1. It 
is crucial to note that light deflection experiments give us privi- 
leged access to y, independently from other post-Newtonian pa- 
rameters. This point is even more important when one thinks 
that cosmological models (Damour & Polyakov 1994, Damour 
et al. 2002) predict deviations of ly- 1 1 of the order of 1 0"*" - 1 0"^ . 

Very long baseline radio interferometry (VLBI) is sensi- 
tive to space-time curvature through the gravitational time 
delay, given by (e.g, Finkelstein et al. 1983) 



GM (\ri\ + ri.k 



(2) 



where r, stands for the position vector of the ith station and 
k the unit vector pointing towards the radio source, both 
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referred to the center of mass of the deflecting body. For 
a typical VLBI baseline between Westford (Massachusetts) 
and WettzeU (Germany) of ~6,000 km, is ~170 nanosec- 
onds (ns) for a source at the Sun's Umb, rapidly decreases to 
~10 ns at 4° away from the Sun, and remains of the order 
of the accuracy of VLBI measurements (nowadays around 
10 ps) even for elongations close to 180° (see Fig.[T)- 
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Fig. 1. Gravitational delay Tg as a function of the elongation an- 
gle ip to the Sun for the baseline Westford- Wettzell. 



VLBI is used on a regular basis since the early 1980's 
for monitoring Earth orientation and estimating station dis- 
placements and extragalactic radio source coordinates at 2 and 
8 GHz. The number of radio sources per session as well 
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as the data recording reliability have drastically improved in 
the last decade. About 4,000 diurnal session files, representing 
more than 5 millions delays, are made available through the 
International VLBI Service for Geodesy and Astrometry (IVS, 
Schluter & Behrend 2007) data base. 

The large amount of data from the permanent geodetic 
VLBI program can provide a number of tests of GR (Soffel 
et al. 1986). In the past years, VLBI data were used in var- 
ious attempts to determine y. Robertson & Carter (1984), 
using less than 4 years of observations, found y consistent 
with GR within 0.005. Robertson et al. (1991), using 10 years 
of observations, estimated a standard error of 0.002. Lebach et 
al. (1995) got 0.9996 + 0.0017 after observations of the relative 
deflection of 3C 273B and 3C 279. Shapiro et al. (2004) obtained 
0.99983 + 0.00026 (statistical standard error) using VLBI obser- 
vations before 1999. The current best estimate of y, however, 
was not obtained with VLBI: it is consistent with GR with an 
error of 2x 10"^ and was obtained by Bertotti et al. (2005) who 
derived it from spacecraft tracking experiments. 

Note that errors reported in the various papers are often 
formal errors obtained from the propagation through the ad- 
justment procedure of an initial SNR-derived standard error 
on the delays. They might therefore not directly compare to 
one another. 

Though the above-referenced works (except Bertotti et 
al.) dealt with deflection of the radio waves by the Sun, it 
must be mentioned that special VLBI sessions were carried 
out to measure the deflection close to Jupiter or other planets 
(Schuh et al. 1988). 

In this work, we estimate y from routine geodetic VLBI ob- 
servations, using the additional 1999-2008 time period with re- 
spect to Shapiro et al. We compare estimates and errors obtained 
over several time spans and using various analysis schemes in 
order to address the accuracy, and to point out some systematics 
and limitations. 

2. Close approaches to the Sun 

A set of 3,937 24-hr geodetic VLBI sessions, consisting of about 
4.5 million delays, will be fully or partly processed in the up- 
coming analyses. During the period that covers 3 August 1979- 
28 August 2008, the VLBI observing schedule included a num- 
ber of radio sources that were observed at less than 15° to the 
sun. As it shows up in Fig. |2] this number was weak before 
1984, quite uniform during 1984-1996. Then it substantially in- 
creased during 1996-2002. It is worth noting that 1992-1999, 
that contains a number of close approaches, is a period of low 
solar activity. Since 2002, the scheduling software at the IVS 
coordinating center was set with a minimal distance to the Sun 
at 15°. Fig. |2] naturally yields several time spans on which the 
analyses can be done: 1979-2008, which is the maximum num- 
ber of available data, 1984-2008, that drops the early VLBI net- 
work, 1996-2002, which shows the highest density of close ap- 
proaches, and 1984-2002, which represents a compromise be- 
tween a high density of close approaches and a large number 
of data. Additionally, we also consider 1979-1999, as done in 
Shapiro et al., in order to check that we are consistent with their 
results. Finally, we would like to address two time spans that 
cover periods of low and high solar activity. It is nevertheless 
difficult to keep the same characteristics (number of sessions, 
number of sources, density of close approaches) for these two 
periods since the VLBI observing program undergoes a con- 
tinuous evolution. We propose the three following time spans: 
1994-1997, 1998-2002a (starts 01/1998, and has approximately 



the same number of sessions and sources than 1994-1997), and 
1998-2002b (starts 07/1998, and has approximately the same 
number of delays as 1994-1997). 

All our VLBI delays have been corrected from delay due to 
the radio wave crossing of dispersive region in the signal prop- 
agation path in a preliminary step that made use of 2 GHz and 
8 GHz recordings. Then, we only use the 8 GHz delays to fit 
the parameters listed in the next section. In the case of targets 
that are close to the Sun, the relevant dispersive regions are the 
Earth's ionosphere and the solar coronal plasma. Although ap- 
proximated, the model for plasma delay correction as a function 
of electronic content and frequency should lead to errors of a few 
picoseconds, following Lebach et al. (1995). (The authors men- 
tioned this magnitude for a period of low solar activity. During 
periods of higher activity, the electronic content can be several 
times higher.) The reader must therefore keep in mind this order 
of magnitude when potential sources of limitation will be listed 
in further sections. Additionally, an error in the solar coronal 
plasma delay correction would lead to a falsified estimate of 7, 
since the plasma-induced deflection would be absorbed there in. 
Rather than a relativistic parameter, y would therefore be simply 
considered as a "deflection" parameter. 
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Fig. 2. Observational history of the sources at less than 30° to 
the Sun. The upper plot gives the Sun spot number (SSN, Clette 
et al. 2007). The right plot displays the deflection angle as pre- 
dicted by GR. 



3. Data analysis and results 

3.1. Global solutions 

We run global solutions over the aforementioned time spans. In 
all these solutions, the Earth orientation parameters and the sta- 
tion coordinates are estimated once per session, y is estimated 
as a global parameter. Source coordinates are also estimated as 
global parameters without global constraint: the sources are al- 
lowed to stay within circles of 10"** rad diameter around a priori 
positions. The choice of the a priori catalogue for source coordi- 
nates is discussed later. 

Now, we quickly go on with some technical characteristics 
of the solutions. The cut-ofF elevation angle is set to 5°. A pri- 
ori zenith delays are determined from local pressure values 
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No. sessions 


No. delays 


No. sources 


Postfit rms delay 
(ps) 


r 




1979- 


-2008 


3,937 


4,386,112 


988 


25.0 


0.99984 + 0.00015 


0.86 


1984- 


-2008 


3,852 


4,348,913 


988 


24.9 


0.99986 + 0.00015 


0.86 


1984- 


-2002 


3,040 


2,857,624 


781 


27.0 


0.99993 ± 0.00017 


0.89 


1979- 


-1999 


2,598 


2,115,509 


723 


27.4 


0.99983 ± 0.00020 


0.91 


1996- 


-2002 


753 


1,024,322 


676 


27.5 


0.99940 ± 0.00022 


0.83 


1994- 


-1997 


650 


849,084 


683 


24.6 


0.99968 + 0.00024 


0.83 


1998- 


-2002a 


650 


953,882 


643 


26.3 


1.00017 ± 0.00032 


0.81 


1998- 


-2002b 


595 


873,827 


616 


26.2 


1.00031 ±0.00035 


0.82 



(Saastamoinen 1972) which are then mapped to the elevation 
of the observation using the Niell mapping function (Niell 
1996). Zenith wet delays are estimated as a continuous piece- 
wise linear function at 20-min interval. Troposphere gradients 
are estimated as 8-hr East and North piece-wise functions at all 
stations except a set of 110 stations having poor observational 
history. Station heights are corrected from atmospheric pressure 
and oceanic tidal loading. The relevant loading quantities are de- 
duced from surface pressure grids from the U. S. NCEP/NCAR 
reanalysis project atmospheric global circulation model (Kalnay 
et al. 1996) and from the GOT00.2 ocean tide model (Ray 1999, 
Petrov & Boy 2004). No-net rotation constraint per session is ap- 
plied to the positions of all stations, excluding HRAS 085 (Fort 
Davis, Texas) and Fairbanks (Alaska) because of strong non- 
linear displacements. (The latter site undergoes post-seismic re- 
laxation effects after a large earthquake on Denali fault in 2003. 
See, e.g., MacMillan & Cohen 2004, Titov & Tregoning 2004, 
2005.) All the calculations use the Calc 10.0/Solve 2006.06.08 
geodetic VLBI analysis software package and are carried out at 
the Paris Observatory IVS Analysis Center (Gontier et al. 2008). 
Results are reported in Table [T] 

Since source coordinates are estimated during the analysis 
process, no significant influence of the a priori catalogue is ex- 
pected on y. To check this, we ran the previous solutions several 
times, using several a priori catalogues. All of them were ob- 
tained after a global inversion of data over 1984-2008, wherein 
the celestial reference frame was maintained by applying a no- 
net rotation constraint on the coordinates on a well-chosen sub- 
set of sources that defines the axes of the International Celestial 
Reference System (ICRS, Feissel & Mignard 1997). Several 
subsets achieve this goal (Ma et al. 1998, Feissel-Vernier 2003, 
Feissel- Vernier et al. 2006, Lambert & Gontier 2009), and en- 
sure an alignment of the output catalogue onto the ICRS within 
0.05 mas. (The latter nevertheless decreases this value below 
0.02 mas.) It finally appeared that the sensitivity of estimated y 
to the chosen set of defining sources and to the a priori catalogue 
is at the level of 10 which is not statistically significant. 

We wondered whether the fit could be improved by remov- 
ing data from sources having a poor observational history (e.g., 
less than 2 observations or observed in less than 3 sessions). We 
therefore ran one more time all the above solutions after having 
downgraded about 200 sources as session parameters and sup- 
pressed the delays from another 100. Final post-fit root mean 
square (rms) and normalized reduced x'' per degree of freedom 
ix^ I f) were not changed significantly. (T\\t I f is output by 
the VLBI analysis software and reflects the goodness of the fit 



of the solution, including all adjusted parameters.) Influence on 
y estimates was only noticed at the level of 10"^, which appears 
to be non statistically significant, provided the standard errors 
reported in Table [1] 

The post-fit rms delay of the solutions ranges 25-28 ps. Such 
a rms corresponds to a rough expected accuracy of 0.27 mas in 
terms of individual source positioning. One can readily see that, 
assuming such a measurement error on the direction of a grazing 
ray one can expect an error 6y not lower than 6yly - 66 jO ^ 
1.5 X 10"^. This is confirmed by the standard errors reported in 
Table □ 

The solutions over 1979-1999, 1979-2008, 1984-2002 and 
1984-2008, that include a large number of sessions and delays, 
and that all have;^'^// larger than 0.86, all result in estimates of 
y consistent with GR within cr = 2 x 10"^. Using the sessions 
after 2002 or before 1979, that do not contain close approaches 
below 15°, makes the estimate of y depart from unity at the level 
of Icr. Incidentally, the solution over 1979-1999 confirms the 
analysis of Shapiro et al. with a slightly lower formal error that 
may originate from a different analysis strategy and a different 
observational data set. 

Although having a similar number of observations or ses- 
sions, solutions over 1998-2002(a,b) bring a standard error sub- 
stantially higher than 1994-1998. Moreover, for 1994-1998, es- 
timate of y appears to be lower than 1, whereas it is larger for 
1998-2002(a,b). It indicates that the bending of sources is higher 
in the latter case. Intense solar activity during the latter period 
could be at the origin of the discrepancy: during periods of high 
activity, the higher electronic content results in a higher deflec- 
tion of radio waves. The absence in the software of a specific 
modeling of solar plasma effects and the strong correlation of an 
uncorrected plasma-induced bending with the relativistic deflec- 
tion prevent from separating these two phenomena. 

3.2. Dependence on the elongation angle 

To address the problem of the elongation angle to the Sun, we 
run several solutions with an increasing cut-off angle, remov- 
ing sources below successive thresholds up to 40°. We apply 
this analysis scheme over several time spans (Fig.|3]l. For 1984- 
2002, a substantial degradation of the estimates occurs beyond 
25°, in agreement with similar tests in Shapiro et al. A bump 
reaches a maximum around 60° and then, estimates of y get 
closer to unity. We run a similar analysis over 2000-2008 be- 
cause it constitutes a data set decorrelated from the one used 
by Shapiro et al. (allowing for the fact that (i) a part of the ob- 
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Fig. 3. Estimates of y for various cut-off of the Sun elongation angle. Horizontal, dashed lines figure ±<t. 
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Fig. 4. (left) Session-wise coordinates of 0229-1-131. (right) Least-squares spectrum. 



served sources and observing antennas are the same in both data 
sets, (ii) the latter contains substantially less sessions than the 
former). The bump also shows up using this data set. 

We also checked what happens at short elongation angles 
over 1984-2002, 1994-1998, and 1998-2002. Below 25°, the 
deviation from unity stays within the error bars with non sta- 
tistically significant variations. For shorter solutions, estimates 
rapidly degrade beyond an elongation cut-off of a few degrees. 
For 1994-1998 and 1998-2002, the degradation occurs in op- 
posite directions. Estimates of 7 appear to be lower than 1 in 
the former case, while they are larger in the latter, consistently 
with the global estimates of y over the same time periods shown 
in Table [T] The possible reason of such differences has already 
been addressed. 

3.3. Approach based on radio source coordinate time series 

Estimating session-wise coordinates of sources can also be a 
mean of looking at a possible deflection when the sources travel 
in the vicinity of the Sun. An uncorrected bending should appear 
as an annual signal in coordinate time series. 

Among the observed sources, only two have close ap- 
proaches below 2° and are observed in more than 500 ses- 
sions. Both cases are similar, and we will only treat the source 



that has the largest observational history: 0229-1-131 (quasar 
4C 13.14). We obtained a coordinate time series using the anal- 
ysis strategy of section [TT] except that y is now fixed to 1, and 
coordinates of 0229+131 are estimated per session. The closest 
approach to the Sun is ~1.5°. At that time, the expected devi- 
ation, following Eq. (1), is ~0.3". Parameter y being fixed to 
unity, this deflection is already corrected and will not show up 
in the coordinate time series. Obtained right ascension and dec- 
lination time series are displayed in Fig. H] The spectrum does 
not show any significant peak at annual period, indicating that no 
extra deflection is detectable. Assuming an hypothetic deviation 
of y - 1 of 2 X 10""*, the incremental deflection would be as drawn 
in Fig.|5] Peaking at ~0.03 mas, it is therefore not detectable in 
the spectrum. It follows that examination of coordinate time se- 
ries for 0229-1-131 can only constrain y to be close to unity at 
approximately the same level of accuracy already obtained from 
global estimates. 

Note that evolution in source structure can show up in coor- 
dinate time series at lower frequencies, as medium or long-term 
patterns (a few months to years), like the slight curvature show- 
ing up in right ascension plotted in Fig.|4j About relations be- 
tween source structure and coordinate time series, the reader 
can refer to, e.g.. Fey et al. 1997, that treats the case of the 
quasar 4C 39.25. 
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4. Discussion and conclusion 

In the above sections, we have used several methods to look for 
radio wave deflection in the vicinity of the Sun, starting from a 
30-yr long routine geodetic VLBI observational data base. We 
interpret this deflection in terms of gravitational bending, as ex- 
pressed in Eq. (1). Using several strategies and various data sets 
covering different time spans, we arrive at the conclusion that j 
is unity within 2 x 10""*. The estimate of y can even reach val- 
ues close to unity by 7 x 10"^ when the time span is limited to 
1984-2002, i.e., to sessions containing observations of sources 
at less than 15° to the Sun. Using longer time spans, although de- 
creasing the formal error due to a larger number of observations, 
makes the estimates depart from unity by about Icr. 

The main factor of limitation is the uncertainty in the deter- 
mination of (global or session-wise) radio source coordinates. 
Causes of this uncertainty have been addressed in various works 
(see, e.g.. Ma et al. 1998, Gontier et al. 2001). The VLBI-derived 
apparent position of a source may change with the global orien- 
tation and shape of the antenna array when the structure of the 
source is extended or not circular. Fey & Chariot (1997), using 
Very Long Baseline Array (VLB A) maps at 2 and 8 GHz, pro- 
vided estimates of the structure delay arising from the extended 
character of the source. In our example of section [331 the struc- 
ture of 0229-1-131 is expected to bring an extra delay below 3 ps, 
let 0.03 mas (see also Ma & Feissel 1997), that partially explains 
the noise level observed in Fig.|4] It turns out that, in absence of a 
direct correction of the delay, based on, e.g., instantaneous maps 
of the source, the accuracy of y estimates from time series 
analysis cannot be better than 10 

Other potential sources of error are the mismodeling of the 
propagation delay through the troposphere, as well as deficien- 
cies in the network (e.g., change of geometry and performances 
from one session to another, dissymmetry between North and 
South hemispheres). The amplitude of the noise that emerges 
from them remains difficult to be precisely quantified at this 
time. It is generally admitted that it is as large as the effect of 
source structure. 

Note also that derivation of radio source coordinate time se- 
ries implies a robust maintenance of the celestial and terrestrial 
reference frames, so that frame effects do not introduce spurious 
perturbations of the estimated coordinates. During the deriva- 
tion process, we checked various analysis strategies and we no- 
ticed that, when the celestial frame is not sufficiently maintained 
(e.g., when too few sources are constrained by the NNR), a semi- 
annual peak could appear at 3cr. In a similar way, fixing the sta- 
tion coordinates to their ITRF values introduces an annual term 
at the same level. These spurious peaks, that could lead to er- 
roneous physical interpretations in the present context, are good 
illustrations of the sensitivity of VLBI to reference frames. 

Already mentioned is the mismodeling of the solar corona 
contribution to light scattering and bending. Although this 
source of error is neglected for geodetic purposes when radio 
sources are observed at large elongations to the Sun, it becomes 
crippling for tests of GR since observers do need to observe as 
close as possible to the Sun. From section 13.11 we tend to con- 
clude that fluctuations in solar coronal plasma limit the accuracy 
of 7 estimates at the same level of the above-listed sources of 
error. 

Thus, various instrumental and modeling errors and analysis 
strategy defects, that cannot be decorrelated and corrected yet, 
explain the current limitation of VLBI for estimating y. 

Compared to the error reported in Shapiro et al., we do not 
consider that we have substantially improved the determination 
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of y. The slight gain in accuracy can be attributed to the ex- 
tra years of data (1999-2008) of which the first 4 years (1999- 
2002) are rich in close approaches, along with the improvement 
of the quality of the VLBI network and observations during this 
time. Our work nevertheless constitutes an independent check, 
and provides some qualitative insights into systematics that show 
up in the analyses of the current geodetic VLBI observational 
database. 

To conclude, we wish to mention that, although current 
VLBI appears to be not competitive with spacecraft systems 
for relativistic experiments, the huge number of VLBI measure- 
ments, in all directions and at a large number of epochs, con- 
stitues an interesting potential for testing other theories than 
the PPN formalism, as for example the scenario of Jaeckel & 
Reynaud (2006) where parameter gamma is replaced by a func- 
tion depending on the elongation angle. 

Acknowledgements. We are grateful to Drs. Anne-Marie Gontier and Peter Wolf 
(Observatoire de Paris) for useful discussions about possible tests. We thank 
Prof. Harald Schuh for his review that helped in improving the paper This study 
could not have been carried out without the work of the International VLBI 
Service for Geodesy and Astrometry (IVS) community that coordinates obser- 
vations, coiTelates and stores geodetic VLBI data. 



References 

Bertotti, B., less, L., & Tortora, P 2003, Nature, 425, 374 

Clette, F., Berghmans, P, Vanlommel, P, et al. 2007, Adv. Sp. Res., 40, 919 

Damour T., & Polyakov, A. M. 1994, Nuclear Phys. B, 423, 532 

Damour, T., Piazza, F., & Veneziano, G. 2002, Phys. Rev. D, 66, 046007 

Feissel, M., & Mignard, F. 1997, A&A, 331, 33 

Feissel- Vernier, M. 2003, A&A, 403, 105 

Feissel-Vemier, M., Ma, C., Gontier, A.-M., & Barache, C. 2006, A&A, 452, 
1107 

Fey, A. L., & Chariot, P 1997, ApJS, 1 11, 95 
Fey, A. L., Eubanks, T. M., & Kingham, K.A. 1997, AJ, 114, 2284 
Finkelstein, A. M., Kreinovich, V. la., & Pandey, S. N. 1983, Ap&SS, 94, 233 
Gontier A.-M., Le Bail, K., Feissel, M., & Eubanks, T. M. A&A, 375, 661 
Gontier A.-M., Lambert, S. B., & Barache, C. 2008, In: D. Behrend and K. 
D. Baver (Eds.), International VLBI Service for Geodesy and Astrometry 
(IVS) 2007 Annual Report, NASA/TP-2008-214162, 224 
Jaeckel, M.-T., & Reynaud, S. 2006, Class. Quantum Grav., 23, 777 
Kalnay, E., Kanamitsu, M., Kisder R., et al. 1996, Bull. Am. Met. Soc, 77, 437 
Lambert, S. B., & Gontier, A.-M. 2009, A&A, 493, 317 

Lebach, D. E., Corey, B. E., Shapiro, 1. 1., et al. 1995, Phys. Rev. Lett., 75, 1439 
Le Poncin-Lafitte, C, & Teyssandier P 2008, Phys. Rev D, 77, 044029 
Ma, C, & Feissel, M. 1997, International Earth Rotation Service (lERS) 
Technical Note 23, Observatoire de Paris 



6 



S. B. Lambert and C. Le Poncin-Lafitte: Determination of y using VLBI 



Ma, C, Arias, E. R, Eubanks, T. M., et al. 1998, AJ, 116, 516 

MacMillan, D. S., & Cohen, S. 2004, In: N. R. Vandenberg and K. D. Baver 

(Eds.), International VLBI Service for Geodesy and Astrometry (IVS) 2004 

General Meeting Proceedings, NASA/CP-2004-212255, 491 
Misner, C. W., Thome, K. S., & Wheeler, J. A. 1973, San Francisco: 

W.H. Freeman and Co. 
Niell, A. E. 1996, J. Geophys. Res., 101(B2), 3227 
Petrov, L., & Boy, J.-P 2004, J. Geophys. Res., 109, 3405 
Ray, R. D. 1999, In: NASA Technical Memorandum, NASA/TM- 1999-209478, 

National Aeronautics and Space Administration (NASA), Goddard Space 

Flight Center, Greenbelt, MD 
Robertson, D. S., & Carter, W. E. 1984, Nature, 310, 572 
Robertson, D. S., Carter, W. E., & Dillinger, W. H. 1991, Nature, 349, 768 
Saastamoinen, J. 1972, In: W. Soren et al. (Eds.), The Use of Artificial SatelUtes 

for Geodesy, Geophysics Monograph Series, Washington, DC: American 

Geophysical Union, 15, 247 
Schluter, W., & Behrend, D. 2007, J. Geod., 81, 479 
Schuh, H., Fellbaum, M., Campbell, J., et al. 1988, Phys. Lett. A, 129, 299 
Shapiro, S. S., Davis, J. L., Lebach, D. E., & Gregory, J. S. 2004, Phys. Rev. 

Lett., 92, 121101 

Soffel, M., Ruder, H., Schneider, M., et al. 1986, In: J. Kovalevsky and V. 

A. Brumberg (Eds.), Proc. lAU Symposium 114, Relativity in Celestial 

Mechanics and Astrometry, Reidel, 277 
Titov, O., & Tregoning, P. 2004, In: N. R. Vandenberg and K. D. Baver 

(Eds.), International VLBI Service for Geodesy and Astrometry (IVS) 2004 

General Meeting Proceedings, NASA/CP-20()4-2 12255, 496 
Titov, O., & Tregoning, P 2005, J. Geod., 79, 196 

Will, C. M. 1993, Theory and Experiment in Gravitational Physics, Cambridge, 
UK: Cambridge University Press 



